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, Abstract. The formal structure of geometrical thermodynamics is reviewed with 

, particular emphasis on the geometry of equilibria submanifolds. On these sub- 

manifolds thermodynamic metrics are defined as the Hessian of thermodynamic 
potentials. Links between geometry and thermodynamics are explored for single 
OO I and multiple component, closed and open systems. For single component, closed 

thermodynamic systems a detailed exposition is given which establishes a clear con- 
' nection between the degeneracy and the scalar curvature of the Weinhold metric 

' (Journal of Chemical Physcis, LV, vol 63) and physical properties such as phase 

o ■ 

' lationship is provided by several specific applications. These include the Ideal gas, 

the van der Waals and the Berthelot gases. For these cases, the degeneracy and 
the scalar curvature of the Weinhold metric are entirely consistent with the actual 



o 



P> I' physical situation. That is, for an Ideal gas no phase transition occurs and there 

■ are no interparticle interactions meanwhile the Weinhold metric is never degenerate 

-4— > ' 

i and has zero scalar curvature. For the van der Waals and Berthelot gases, both 

display a phase transition and experience interparticle interactions, the Weinhold 
metric is generally degenerate along a sub-manifold of co-dimension one and has 
, non-zero scalar curvature. For multi-component closed and open systems the Gibbs 

free energy is employed as the thermodynamic potential to investigate the con- 
nection between geometry and thermodynamics. The Gibbs free energy is chosen 
for the analysis of multicomponent systems and, in particular, chemical reactions. 
Major emphasis is focused on a detailed examination of single chemical reactions 
in a multi-component closed system. Then the approach is extended to considera- 
tion of / independent chemical reactions. For single chemical reactions, the general 
Gibbs metric for the Ideal gas mixture is provided. Specific applications include the 
isothermal and isothermal-isobaric cases for a simple synthesis and single displace- 
ment chemical reactions. For these simple systems, results suggest an intriguing 
relationship between non-ideal interparticle interactions and phase transitions. Fi- 
nally, the Gibbs metric is provided for multi-component open systems, including 
both ideal and non- ideal solutions. 



1 



2 



1. Introduction 



In the traditional approach to present the basic structure of homogeneous ther- 
modynamics, it is customary to fix a set of variables describing the state of a system, 
the processes going on in the system and interactions with the outside world. Such 
a set usually includes the internal energy U of the thermodynamic system. A set of 
values of these functions form the extended state space of the system (referred to 
as the energy-phase space in [8], and thermodynamic phase space in [11]). All but 
one of these variables, named thermodynamic potential in this representation, are 
set into couples {y,x) of intensive and extensive variables in such a way that to each 
extensive variable x there corresponds an intensive variable, y, and the infinitesimal 
work (change of energy U or a chosen thermodynamic potential) related to the change 
in the extensive variable x is. 



These couples are often collected into larger groups, corresponding to the tensorial 
type of the process they describe or to the process in whose description they partici- 
pate. Collecting all such pairs, the first law of thermodynamics in its geometrical 
form postulates that during the process, the change in the internal energy, U, is given 
by integration over the trajectory in the state space of the one form. 



Here dQ is the heat change one form. The second law of thermodynamics in the 
formulation of C. Caratheodory [6,7] states that the form. 



has an integrating factor (see also [3]). After some effort [6,7], this integrating factor 



dW — ydx 





is determined to be ^ 



(T is the absolute temperature). Thus, 



dQ = TdS 
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where the new state variable S is the entropy. 

The couple of variables {S,T), (entropy/temperature), plays a special role in this 
formulation. If one lists all the extensive variables Xi (including entropy) and the 
corresponding intensive variables yi, the infinitesimal change of the internal energy of 
the system is given by, 

k 

dU = yidx" 
1=1 

Thus, if the thermodynamic phase space in variables {U;x\yi) is denoted as P, 
there is a one-form 9 defined by the choice of the process and the variables related to 
it, namely, 

k 

e^dU-^Vidx' (1.1) 

i=l 

Processes that might occur in the system should be such that, along the curve 
t — > r{t) — {U{t),x^{t),yi{t)), 9{r'{t)) — 0. Thus, in this geometrical situation the 1- 
form 9 defines the contact structure on P and all the physically admissible processes 
should be integrable curves of the contact distribution Dp — ker{9p), with p e P, of 
this structure [6,7]. 

1.1. Geometrical thermodynamics. Geometrical interpretations of equilibrium 
thermodynamics have proved that state space is endowed with a canonical contact 
structure that underlines the first law of thermodynamics. Different representations 
of this structure in a canonical D'Arbois-chcirt are related to different forms of 
the law of conservation of energy which can be expressed through the internal energy, 
entropy, Helmholtz free energy, or other extensive variables [11,15]. Hermann [10] 
and later Mrugala [11] argued that expended phase space of a homogeneous thermo- 
dynamic system, endowed with the contact structure, is the natural geometric space 
for descriptions of equilibrium thermodynamics. Until now, explicit application of 
this geometric analysis to gain insight into, among other things, the critical behavior 
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of real chemical systems has not been presented. Applications of this geometric ap- 
proach to the analysis of simple thermodynamic systems is the focus of the present 
study. 

To begin, R. Hermann [10] and later R.Mrugala [11] defined the extended state space 
of a homogeneous thermodynamic system as a (2k+l)-dimensional manifold P 
endowed with the contact structure given by a differential 1-form 6 such that 

^A(d^)VO (1.2) 



where 9 is called the contact form. This condition is equivalent to the property of the 
smooth subbundle D C T{P) being as far from integrable as possible [1]. 
Locally, the association p — > ker{9p) defines a 2k-dimensional distribution 



independent of the choice of 9. 

Moreover, it is possible to show that replacement of 9 hy f{p)9, with some function 
f{p) positive at all points of domain of 9, does not violate this condition. Any contact 
form 9, in an appropriate local canonical chart C of variables x'',yi,i — 1, . . .k), 
called the D'Arbois chart, is expressed by [1], 

fc 

9^dx° -J^Vid^'- (1-3) 
1=1 

In such a canonical chart, the distribution D at the point p is generated by the 
vector fields, 

d d d ,^ ^, 



and the differential of the form 9, d9, defines on each hyperplane. Dp, the symplectic 
structure {Dp,d9). In local coordinates d9 has the canonical form. 
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k 

d9p = dx"^ A dyi 
1=1 

Replacement of 9 by f{p)9 leads to the replacement of d9 by df A9 + fd9 which, after 
restriction to a hyperplane Dp, becomes fd9. Thus, contact structure alone defines 
only conformally symplectic structure on the distribution D. 

On the manifold P there exists a unique smooth vector field Y called the Reeb 
vector field [5] such that, 

9(Y) = 1 iY(d9) = (1.5) 

In particular one gets the canonical splitting, 

Tp{P) ^ Dp®keTd9p (1.6) 

of the tangent bundle of P into the direct sum of two subbundles, the first being 
the subbundle of horizontal vectors of the distribution D, while the second being the 
characteristic subbundle of the form d9 [5]. Correspondingly, the cotangent bundle 
T*{P) splits as well. 

In the canonical D'Arbois chart, the Reeb vector field is just. 



Couples of variables {x\ Hi) denote pairs of independent parameters (x*) and corre- 
sponding conjugate variables (y^) with respect to the chosen thermodynamic potential 
x°. Examples of {x\yi) pairs are: (1) temperature and entropy {T,S); (2) pressure 
and volume {p,V); (3) mole number of i-th component and corresponding chemical 
potential {Ni, extent of reaction and corresponding affinity {i,A). 

1.2. Thermodynamic Equilibrium. Thermodynamic equilibrium is a key notion 
in thermodynamics. In particular, 

in all systems there is a tendency to evolve toward states in which the properties 
are determined by intrinsic factors and not by previously applied external influences. 
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Such simple terminal states are, by definition, time independent. They are called 
equilibrium states [4]. 

In a state of thermodynamic equilibrium, intensive variables are functions of the 
extensive variables, namely 

Vi = Viix^) (1.7) 
Then, choosing the internal energy U as the potential, the form 9 becomes 

k 

9^dU-^yi{x^)dx' (1.8) 

i=l 

Along the contact distribution D — ker{9), 

dU ^^yi{x^)dx' (1.9) 

i 

Thus, the relation yi{x^) — ^ exists just on the maximal integral submanifold of 
the contact manifold P. 

Denoting $ as a generic thermodynamic potential, it is known that equilibrium 
states belong to a maximal integrable surface of contact form 9 in the space P deter- 
mined by the choice of k independent variables x^ and by the thermodynamic potential 
$(x*) as a function of these variables (constitutive relations). Another choice of in- 
dependent variables together with some other specification of $(a;*) leads to another 
equilibrium surface corresponding in general to another constitutive relation. The 
core of the present study focuses on Legendre submanifolds defined to be maxi- 
mal integral k-dimensional submanifolds of P on which the Pfaff equation ^ = holds 
[1,11]. The standard approach to locally defining such a submanifold, <S$, in terms of 
a generating function, $, is given by the following theorem: 



Theorem 1. (V.Amold, [1]). 
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For any partition lUJ of the set of indices (l,...,k) into two disjoint subsets I.J and 
for a function ^{yjjX"^) of k variables yi with i E I and with j G J, the following 
equations, 

define a Legendre submanifold 5$ of a contact manifold (p2fc+i^ 

Conversely, every Legendre submanifold of {P'^^^^^D), in a neighborhood of any 
point, is defined by these equations for at least one of 2^ possible choices of the subset 
I. 

In the special case in which $ is a function of only the independent variables 
the Legendre submanifold (submanifold of equilibria states) <S$ is given 

by, 

= {{^,x\y,,t = 1, . . . A;) G P|$ = = e I) (1.11) 

On the integral submanifold 5$, the function $(a;*) can be defined in terms of the 
variables x\ yi as, 

^{x\...,x') = Y,^^yi{x\...,x') (1.12) 

i 

In most cases, (1.12) is homogeneous of order one and satisfies the Euler equation, 
i.e. 

$(Ax\ Ax^') = \^{x\ ...x'') (1.13) 

(see [9]). The expression in (1.12) leads directly to the Gibbs-Duhem relation between 
variables along the integral submanifold S^. Indeed, taking the differential of $, we 
obtain 



Thus, the basic contact condition ^ = is equivalent to the Gibbs-Duhem consti- 
tutive relation, i.e. 
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x'dyi = 



In order to apply the contact formalism to equilibrium thermodynamics, all vari- 
ables Hi, x\ i — l,...,k, must be identified with the thermodynamic parameters 
such that 9 — dx^ — yidx^ satisfies the first law of thermodynamics. These parame- 
ters only have physical meaning on the k- dimensional Legendre submanifold defined 
by the Pfaff equation ^ = 0. In this context, = $ is a thermodynamic potential 
function of the independent variables x^,...x''; and yi,---,yk are the corresponding con- 
jugate parameters with respect to the potential. On the Legendre submanifold, the 
parameters ?/i = [13]. 

One of the goals of the present study is, in the context of geometrical thermody- 
namics, to show how the choice of extensive variables and thermodynamic potential 
is crucial for a reasonable physical interpretation of geometrical objects. 

In the present development, two sets of variables are chosen. In the case of a single 
component system. A; = 3 and is identified with the internal energy, while, x^, 
x"^, x^ are equated to the independent variables S, V, N, respectively. Likewise Pi, P2, 
P3 are the corresponding conjugate variables T, -p, ^, respectively, and the contact 
form becomes 6 = dU — TdS + pdV — ixdN. This is an extension of previous work 
[20] with the addition of a more detailed exposition on the relation between geometry 
and thermodynamics. 

In the case of an r-component system, k = r + 2 and x^ is identified with the 
Gibbs free energy; x^, x^, x^^.-.^x^^"^ with the independent variables T, p, Ni,...,Nr; 
and pi, p2, p3,...,Pr+2 with the conjugate variables -S, V, /Xi,...,//^, respectively. Then, 
the contact form becomes 9 — dG + SdT — Vdp — jXidNi — ... — jirdNj.. 

1.3. Thermodynamic Metric. A thermodynamic metric ?7$ defined by the consti- 
tutive relation $ = <l>(x*) on the Legendre submanifold of the contact structure 9 
has the form( [14,20]), 



For the case where $ is the internal energy, U, the metric, rfu, is called the Weinhold 
metric([21]). When $ is the entropy, S, the metric, 775, is called the Ruppeiner 




(1.14) 
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metric([r7]). Here a new metric, r/c, is introduced for the case where $ is the Gibbs 
free energy, G. 

Remcirk 1. A thermodynamic metric, 77$, of the form (1.14) is induced on <S$ by the 
following symmetrical tensor ([20]), 



Up to a conformal factor, this tensor is the symmetrical tensor in P annihilating the 
Reeb vector field, Y, of structure 6, (Y = gfo), and is invariant under substitution of 
indices, i. fj is obtained as the sum of symmetrical tensors in the 2-dimensional sub- 
spaces D* of D* spanned by pairs of covectors {dx\ di/i) of thermodynamic conjugate 
variables [20]. 

Moreover, thermodynamic metrics are generally degenerate and non-definite. In 
particular the Legendre submanifold or (as it is referred to hereafter) the thermo- 
dynamic state space, is the union of domains where these metrics have different 
signatures separated by the submanifold (generically of codimension one) of states 
where these metrics are degenerate [20] . 

1.4. Scalar Curvature of the Thermodynamic Metric. Now consider a ther- 
modynamic potential, $ = $(x*), a function of the independent variables, x*, and 
calculate the scalar curvature of the corresponding thermodynamic metric defined on 
the Legendre submanifold 5$. According to expression (1.14), 




(1.15) 



(1.16) 



Christoffel symbols for this metric are given by (see also [18]), 




(1.17) 



m 



where rjij^rn = 
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It can be shown that the curvature tensor of the metric rjq, is given by, 

"^ij-fe — ^ki,j ~ ^''ji,k + Tj-pr^j — r^pF^j = -^{Vij,mVsn,k — Vsn,jVki,m)'n^^V^^ (1-18) 

Therefore, the Ricci Tensor of metric rjq, is given by, 

T^ik ~ '^ijk ~ '^iVij,rnVsn,k Vsn,jVki,m)V ) (-'-•-'-^) 

and the scalar curvature TZr,^ by (see also [18]), 

= '^ikV''' = ^{Vij,mVsn,k - Vsn,jVki,m)v'^Vv''' ■ (1-20) 

Remark 2. Consider the scalar curvature of the metric defined on two-dimensional 
integral surfaces <S$. Components of the Ricci tensor are given by([20]), 

■^/// \2 \11 22 / \ 12 22 

= ^(((?72i,J -Vn„V2iJV V + (r/2i,i?72i,2 - ^11,1 ?722, Jr/ r] + 

((^2iJ^-mi„r/22j(ry^^)^) (1.21) 

7^12 = 1^21 = ^{{VU„V21,^ - {V2l„fW^V^^ + (^711,1^22,2 - r/21,2?7ll,J(^^^)^ + 

{vu,,m2,, - imiJ'WW') (1-22) 

T^22 = ^((?721,J^ - V22„VuJ{V^^f + {V2l„VU,2 ' ?7ll,i^22, J?7^^??^^ + 

((r/2ij'-r/i2,r722j^''r7'') (1-23) 
It is straight forward to show that, 

7^ll?7'' = ^22?/'' (1.24) 
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7^ = 2(7^nr/" + 7^l2r]^') (1.25) 

These general expressions for curvature are employed to characterize the thermody- 
namics of several systems. 

1.5. Closed and Open Thermodynamic Systems. Consider the thermodynamic 
state of a system as a function of a certain number of independent variables such as 
the entropy, S, volume, V, and number of moles, Ni,...,Nr, of r components. Any 
function expressed in terms of these variables is a state function of the system. The 
First Law of Thermodynamics or Conservation of Energy postulates the existence 
of a state function, called the energy function, such that the change in internal energy 
of the universe, given as the sum of the energies of our system and of the surroundings, 
is always constant, namely([16]), 

= dU^Yiiv. ~ dUgyg^ -\- dUsurr. (-^•^^) 

This expression states that, 

dUgys. — dUgurr. (^■^'^) 

The minus sign indicates a loss of energy by the surroundings. Denoting dsU — 
—dUgurr. as the change in energy supplied to the system by the surroundings. Con- 
servation of Energy dictates, 

dUgyg, ^ dU ^ dsU or djU^O (1.28) 
Subscript / indicates the energy change of the system. Note, djU — 0, is equivalent 

to dUuniv. = 0. 

The Second Law of Thermodynamics or the principle of Entropy Production 
postulates the existence of a state function, called the entropy function, which pos- 
sesses the following properties: the entropy is an extensive variable, and the change 
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in entropy dS can be separated into the flow of entropy, (IeS. due to interactions witli 
surroundings and a term, diS, corresponding to entropy cliangcs in tlie system ([16]). 
Tliat is, 

dS = dES + diS (1.29) 



where djS is denoted as the entropy production. djS is always non-negative, zero for 
reversible processes and positive for irreversible ones. 

For closed systems, conservation of energy in (1.28) can be expressed as, 

dU ^dQ + dWM ^dQ -pdV (1.30) 



with pressure, p, normal to the surface. For open systems, 

d?/ = + dWM = - pdV (1.31) 



where d'if is the infinitesimal rate of change of heat transfer and exchange of matter 
of the system with the external environment [16]. These expressions are employed 
to examine the geometries of three different thermodynamic systems. First, work 
done on single component thermodynamic systems is reviewed while stressing the 
importance of choosing the right thermodynamic variables suitable for a particular 
situation. For a one component thermodynamic system physical interpretations are 
deduced from geometrical objects such as the degeneracy and scalar curvature of the 
Weinhold metric. Choosing the Gibbs free energy as the preferable thermodynamic 
potential certain physical aspects of chemical behaviour can be described through the 
geometry. This approach is also used in studying chemical reactions in multicompo- 
nent systems. 

1.6. System 1: Single Component Closed System. Consider a closed system 
containing a single component in the absence of an external field ([20]). The energy 
supplied by the surroundings is derived from the sum of the heat transfer, dQ, and 
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mechanical work, cIWm- In this case the entropy production diS = 0, and the entropy 
of the system is given by, 

dS^dES^^ (1.32) 



Therefore, in molar form the 1-form of the energy, du, is given by, 

du = dQ + dWM = Tds - pdv (1.33) 

1.7. System 2: Multi-Component Closed System. Consider a multi-component 
system in which changes in internal energy can occur due to chemical reactions. In 
this case, the entropy production djS is given by [16], 

djS = pi 

where A is the affinity of the chemical reaction related to the chemical potentials jii 
hy A — —^j^jJiiVi- The Vi are the stoichiometric coefficients and ^ is the extent of 
reaction. Note, for a single chemical reaction the entropy of the system is given by, 

dS = dES + diS=^ + pi (1.34) 
Now introduce another state function, the Gibbs free energy, G, defined by, 

G^U -TS + pV 
In terms of this function, conservation of energy can be written as [16], 



dG^dQ- TdS - SdT + Vdp = 



-SdT + Vdp - Adi 



(1.35) 
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While this expression is central to the present study, useful geometrical tools are 
also introduced allowing a more general treatment of the case of I independent chem- 
ical reactions. In this context, (1-34) and (1.35) can be restated as, 

dS = ^ + J2j^d^n, (1.36) 

n 

and 

dG = -SdT + Vdp-J2 ^nd^n (1-37) 

n 

where n — 1, is the number of chemical reactions that occur. 



1.8. System 3: Open Systems. Consider an open system in the absence of an 
external field without entropy production. The expression in (1.29) can be generalized 
to consider changes in the number of moles Ni,...,Nr [16], viz. 



dS = dsS = 




(1.38) 



i 



where d'if — TdS + J2i l^idNi is the energy flow due to heat transfer and exchange of 
matter. The corresponding Gibbs free energy is given by, 

dG = -SdT + Vdp + Y^ l^idNi (1.39) 

i 

with i — 1, r, the number of moles of each reaction component. 

2. System 1: Single component system 

Consider a 7-dimensional thermodynamic phase space P of variables 
([/, {S, T), (V, —p), {N, ji)) with the contact 1-form given by. 



e^dU - TdS + pdV - iidN 



(2.1) 
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where is the number of moles of the component and fi its chemical potential. Next, 
consider a 3-dimensional Legendre submanifold Su{S, V, N) of this system defined by 
the constitutive relation, 



U = U{S, V, N) (2.2) 

By homogeneity of degree one of the internal energy, consider the molar form of the 
constitutive relation (2.2) and obtain the following constitutive relation for a closed 
system with a single component [4] , 

u = u{s, v) (2.3) 



The differential is given by (1.33), i.e. 

du — Tds — pdv 



Introduce the following thermodynamic parameters: 
(1) Cy is the heat capacity at constant volume: 



Cv = T{—)y, (2.4) 



(2) Cp is the heat capacity at constant pressure: 



C, = r(— )„ (2.5) 



(3) a is the thermal coefficient of expansion: 



1 dV 

(4) is the isothermal compressibility: 



a = (2.6) 
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For the molar case the above parameters are represented in lower case type. The 
following expression relates the above molar parameters: 

Cp- c^ = (2.8) 

Considering the expression in (1.14), the Weinhold metric rju is defined on the 
two-dimensional integral surface Su by [12], 

'"• dx'dxi c„ l-Ji! ^] ^ ' 

\ Kt VKt / 

It follows, if T 7^ 0, the Weinhold metric is degenerate along the curve 7^ given by, 

(|). = (2.10) 



which is presented in one of two forms: p — p{v) or T — T{y) 
The main points of focus are as follows [20] : 

1) The equilibrium surface is the union of regions where the Weinhold metric has 
different signature separated by the curve 7^ where the metric is degenerate; 

2) The critical point (pc, T^i Vc) of the system is the extremum of the functions 
p = p{v) and r = T(i;); 

3) Along the curve of degeneracy, 7^, a first order phase transition seems to occur; 

4) Scalar curvature of the Weinhold metric is strongly influenced by parameters 
related to non-ideal inter-particle interactions within the system. 

Remctrk 3. For a two dimensional state space, the determinant det{r]^) and the 
scalar curvature TZ oi a thermodynamic metric 77$ are inversely related [20] , 



TZ — —— — z r:;det 

4:det[r]^y 



f VU VU,! ^11,2\ 

^12 ?7l2,l r]l2,2 
\rj22 ?722,1 ?722,2/ 



(2.11) 



In point 4) above, it was noted that the scalar curvature of the Weinhold metric is 
strongly related to interparticle interactions in the system. Moreover, as the deter- 
minant of the matrix approaches zero, point 3) implies that the system approaches a 
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phase transition at which point the scalar curvature TZ of the metric goes to infinity. 
Thus, the connection of points 3) and 4) suggests the following interpretation: 

5) If the system approaches a state "close enough" to the curve of degeneracy, 7^, 
the scalar curvature of the metric goes to infinity. Physically this is consistent with 
a relevant increase in inter-particle interactions between the reactant and product 
species when the system approaches a phase transition. 

Such behavior suggests an intriguing relationship between degeneracy, scalar cur- 
vature and inter-particle interactions. In particular, this suggests a geometrical condi- 
tion for a phase transition might be degeneracy (or infinite curvature) of the Weinhold 
metric r]ij^. The following examples support this suggestion. 

Example I: Ideal gas. 

Given the equation of state for an Ideal gas, pv — RT, the Weinhold metric, ry^^, is 
given by [19], 

In this case Cp — = R, and [20] , 

V vT RT'^ 

detiVijJ ^^,{Tcp-pv)^^^^>0 (2.13) 

The above metric is positive definite on the costitutive surface <S„. Thus, for an 
Ideal gas, the energy metric is never degenerate except for trivial cases, i.e. T = 0. 
This lack of degeneracy is consistent with the characteristics of an Ideal Gas (i.e. it 
does not display a critical point and therefore does not exhibit a phase transition). 
Moreover scalar curvature of the Weinhold metric is zero, [17,20], i.e. 



= (2.14) 

Ruppeiner [17] was the first to suggest that zero curvature might be evidence for 
the absence of inter-particle interactions. Their absence is precisely the case for an 
ideal gas. Table 1 gives parallels between the geometric and thermodynamic features 
of an Ideal Gas. 
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Geometry 


Thermodynamics 


No curve of degeneracy 


No phase transition - no critical point 


Zero scalar curvature 


No inter-particle interaction 



Table 1. Ideal Gas: Geometry- Thermodynamics 



Example II: van der Waals gas 

The equation of state for the van der Waals gas is given by, 



{p+^){v-b) = RT (2.15) 



y2 ' 



where a and b are positive constants. Expression (2.15) provides a more realistic 
representation of the actual behavior of real (non-ideal) gases by introducing the 
additional positive constants a and b, characteristic of the particular gas under con- 
sideration. The factor {v — b) indicates the excluded volume of the molecules, while 
the factor is the "interaction" term. 

For the van der Waals gas, the Weinhold metric is given by [20], 

T TR 



Viju - I TR / TR _ 2a\ I ^"^'^^^ 



which is degenerate along the curve 7^ written in the following forms: 



s = s{v) = c,[(2 + ^) In {v-b)+ In (^; 



p{v) = {v-2b)-, T{v) = ^^^3 



Note, in the hmit a = 6 = 0, the ideal case is recovered with the metric in (2.12). 
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Figure 1. p-T phase boundary 

Taking the derivatives of the last two expressions p — p{v) and T — T{v) and setting 
them to zero, the critical point is obtained, 

Moreover, for the van der Waals gas the scalar curvature of the Weinhold metric 
'^Viju is given by [20], 



T^n = —, 5 TTT (2.18) 



c„(pt;^ — af + 2a6)^ 

In general, the scalar curvature 7^r;ij„ as a ^ and as the system approaches 
the degeneracy curve, 7^,,^^^ — > oo. On the other hand, expression (2.18) does not 
vanish if the parameter b ^ 0. So while the scalar curvature of the Weinhold metric 
is not related to the excluded volume, it is strongly influenced by the parameter a, 
that includes non-ideal interactions in the system. This finding is entirely consistent 
with the physical behavior of the van der Waals gas which exhibits both a critical 
point and phase transition. Indeed in the (p — T) plane, the following solutions are 
obtained [20], 

,-1 3 



One of the these solutions is the Pressure- Temperature Phase Boundeiry (Fig.l). 
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Table 2 summarizes parallels between geometric and thermodynamic features of 
the Van der Waals gas. 



Geometry 


Thermodynamics 


Curve of degeneracy 7,, 


phase transition (see Fig.l) 


Extrcmum of 7^ 


Critical point 


Non-zero scalar curvature 


Inter-particle interaction 



Table 2. van der Waals: Geometry-Thermodynamics 



Example III: Berthelot gas 

As a third example consider the Berthelot gas with the equation of state, 

{p+^){v-h) = RT (2.19) 

where a and b arc positive constants (analogous to the van der Waals gas). 
The Weinhold metric for the Berthelot gas is given by [20], 



T_ 1 t RT _| a_\ 

„.. _ I c„ cv^{v-b) ^ Tv^> 1 (oor\\ 

Vi3u-\ ±(^J_:_g_\ RT^v''-2a{v-bf , T( R , a ^2 t 

whose degeneracy is [20], 

RT^v^ - 2a{v - bf = 2p^v%v - hf - aR{v - 2hf = (2.21) 
In analogy to what was done for the van der Waals gas, it follows that 

, „ , , , , aR , 1 , 80, . 1 . 
(.„p„TJ = (36,±(^)..±(— ).) 

and that the scalar curvature of the Weinhold metric, TZr^j^ , is given by [20] , 

^ ^ {T%''Rc,L{c,,v) +T^v'RaQ{c,,v) + a^Wjcv)) 

'^viju « clT^v{RT^v^ - 2a{v - byy ' ^ ' ^ 



where 
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L{cv, v) — (2c„ — R)v^ — 2>Cvbv + c„6^ 
Q(c^, v) = -Rv^ + mbu^ - SRb'^v^ + {Rb^ + c„ + R)v^ - b{b - 2v){R + c^) 

and 

W{cy, v) = -Rv'^+ARbv^-6Rb'^v^+{2cy+R+ARb^)v^-{Scy+3R+Rb^)bv^+{12cy+3R)b'^v^- 

-{8cy + R)b^v + 2cyb^ 

The scalar curvature, TZruj^ , goes to zero as a — 0, and is strongly influenced by this 
parameter which corresponds to non-ideal interactions in the system. Furthermore, 
TZriij^ —>■ oo as the system approaches a phase transition. Once again, degeneracy of 
the Weinhold metric and non-zero scalar curvature are consistent with the character- 
istic physical behavior of the Berthelot Gas. Parallels displayed in Table 2 for the 
van der Waals gas are equally applicable to the Berthelot gas. 

3. System 2: chemical reactions in a closed system 

3.1. Single chemical reaction. Consider a closed system comprised of r compo- 
nents among which chemical reactions can occur. First, we focus our attention on 
single chemical reactions and then introduce multi-component reactions. In a closed 
system, any change in the masses of the components will occur only from a chemi- 
cal reaction. Thus, denoting the mass of component i by rrii, with i — 1, ■■■,r, the 
infinitesimal change in mass can be written as [16], 



(3.1) 
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where Mi is the molar mass of component i. The principle of conservation of mass 
for a closed system is expressed as [16], 

dm^Y^ UiMid^ = (3.2) 

i 

with m = 'Ylii'^i- The equation '^■ViMi = is referred to as the stoichiometric 
equation. 

Alternatively, rather than the component masses it is more convenient to consider 



dm 

tesimal change in the mole number of the i component, can be expressed as 



the number of moles Ni,...,Nr involved in the reaction. Since = dNi, the infini- 



dNi — Uid^ 

Let be the number of moles of component i in the initial state of the system. When 
a reaction occurs, as indicated by the stoichiometric coefficients Ui, the variations of 
the number of moles of each component Ni are not independent. This can be expressed 
as [9], 



= ... = = ... = (3.3) 



where the extent of the reaction, ^, is an extensive variable just like the number 
of moles. Integrating and taking ^ = as the initial state of the system, we obtain 
[9], 

Ni = N^ + Ui^ i = l,...,r 



In this context, the Legendre submanifold Sq can be defined by the constitutive 
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relation, G = G{T,p, ^). Restriction of the contact 1-form, 6 = dG+SdT — Vdp+Ad^ 
to the submanifold Sq provides, 

dG = -SdT + Vdp - Ad^ (3.4) 

Thus, the general metric rjij^ — g^ig^j , where and are the extensive variables, is 
given by. 



Vija ~ 



"t 

aV -krV ArV 



( aV -ArS \ 



(3.5) 



\-A,S A^V -mrj 



where —A = Aj.G is the Gibbs free energy of the reaction and Aj.S and Aj.V are 
the entropy and volume of reaction, respectively. Here, the affinity and Gibbs free 
energy of reaction arc used intcrchangebly. 

Naturally, as a reaction takes place, the chemical potential of the components 
varies and so does the affinity of the reaction. At constant temperature and pressure, 
the system is at equilibrium whenever the affinity A = 0. Since (^)r,p < 0, the 
determinant of the matrix in Eqn.(3.5) is given by, 

detv^Jo - -^^i^h,, + Y^A^Vr + VkriA^S - ^A^V]' > (3.6) 

Of primary interest is what type of information is provided by the degeneracy and, 
in some simple cases, by the scalar curvature of the metric ryjjg. As an example 
consider the three-dimensional case of the Ideal gas mixture. 

Example I: Ideal gas mixture. 

Consider a simple reaction in which substance A converts into substance B. Starting 
with one mole of A, the relative amounts at some later point in the reaction are 
Na = 1 — ^ and Nb — ^ with ^ e [0, 1]. Therefore, the Gibbs free energy can be 
written in terms of the extent of reaction as. 



G = (1 - O/^A + ^f^B 
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Now, the chemical potential of an ideal gas mixture is given by /Xj = /if(T) + 
i?T In(p-) where i = A, B and the superscript 6 indicates some standard state at 
pressure p^. This reaction resulting in conversion of the ideal component A into the 
ideal component B is analogous to a homogeneous mixture of two ideal components, 
in that the two components are mixed but do not interact. 

Considering that = (1 — and ps = ip, where p is the total pressure, the 
Gibbs free energy can be written as, 

G={1- OA + ^I^B + RTH^) + RT[{1 - ln(l - + ^In^] 

where -^^[(l — ^) ln(l — + Cli^C] = ^Gmix is the Gibbs free-energy of mixing. 
Thus, the metric (3.5) becomes. 







\ dT dT ^ -"-^"li-^J ^ 



p 



RT 

n2 



RT 



/ 



(3.7) 



where —ArS 



^ + i?ln(^), ArV = and (|f)r,p = -J^- The determi- 
nant in expression (3.6) reduces to. 



dT 



det,«„ = ^H^iA - A) + RHj^^)? - + (1 - 0# + #1) 

(3,8) 

If the chemical potentials of the two ideal components in the standard state are 
explicitly known, useful general information could be extrapolated from the Gibbs 
metric, its degeneracy and its scalar curvature. Since in general this is not the case, 
our analysis is restricted to the 2-dimensional isothermal case and to the 1-dimensional 
isothermal-isobaric case. 

3.2 Isothermal single chemical reaction. 

When the temperature is kept constant during the chemical reaction, dT — and 
the expression in Eqn. (3.4) reduces to. 



dG = Vdp - Ad$, 



(3.9) 
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Thus, the metric (3.5) reduces to, 



-krV A,.V 



with 

,dA 



detrj,j,^kTV{^),-{ArV)' (3.11) 
In the case of an ideal gas mixture, the metric of Eqn.(3.10) becomes. 



_ju_] (3-12) 

with determinant det rjij^ — —^^z^ which is always different than zero (except for 
trivial values of some thermodynamic parameters) . This implies that the Gibbs metric 
is, in general, never degenerate for an Ideal mixture. Thus, the physical interpretation 
of this result is that there is no critical behavior displayed by an Ideal Gas Mixture. 
Moreover, the scalar curvature of the metric (3.12) is zero. Indeed, since the inverse 
of riij^ is given by, 

€ - ("o* 4il) (3.13) 

the third derivatives of the Gibbs potential are given by, 

2,RT 

Vu„ = Vu„ = Vi2„ = V2i„ = (3.14) 

i?r(2g - 1) 

V22„ = 7712,2 = V21„ = 0, 7722,2 = ^2(1 _ (^'l^) 

Therefore, the components of the Ricci tensor TZij, (1.21) to (1.23), for an isothermal 
ideal mixture are all zero, namely 
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Using the expression in Eqn.(1.25), 

= (3.16) 

Obviously, this result is consistent with the fact that the two Ideal components 
when mixed do not interact, and it is essentialy consistent with the features of the 
single-component Ideal case. This strongly suggests that even in the context of chem- 
ical reactions in closed systems, non-zero scalar curvature might provide useful infor- 
mation regarding interactions between components. Although beyond the scope of 
the present study, this is an interesting path to pursue and, due to the complexity of 
the system, will require the use of numerical mathematics. 

3.3 Isothermal-isobaric single chemical reaction 

It is interesting to note that in the case of constant temperature and pressure, the 
change in Gibbs free energy is given by, 

dG = -AdC (3.17) 

In this one-dimensional case, important information can be gleaned from examination 
of the convexity of the Gibbs free energy function. Consider the condition. 



d'G dArG 



de d^ 

For an ideal gas mixture [9], 

d^G RT 



= (3.18) 



which implies that the Gibbs free energy is a convex function of the extent of re- 
action. An example of such a function is displayed in Fig. 2, (see [22]), and defines 
the condition of stability. Initially, the Gibbs free energy decreases. As a reaction 
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^ Extent of reaction, ^ 

Figure 2. Ideal Gas Mixture: G vs. ^ 

proceeds, the Gibbs free energy of the system continues to decrease until it reaches 
a minimum value. At equilibrium (constant T and p), the Gibbs energy is at the 
minimum, and A^G is equal to zero. At greater extents of reaction the Gibbs free 
energy is greater than zero and increases. 

This well-known result suggests for a non-ideal mixture the change in sign and 
vanishing of the second derivative of the Gibbs function might provide some insight 
into the stability of the system. In particular, consider the case in which the Gibbs free 
energy is not a simple convex function of the extent of the reaction everywhere, but 
rather, first convex, then concave, then convex again, as shown in Fig 3. Points on the 
curve where changes between convex and concave behavior occur are inflection points 
which satisfy the expression in Eqn.(3.18). Naturally, a shift of the equihbrium state 
from one local minimum to another constitutes a first order phase transition induced 
by a change in the extent of reaction [4]. For a chemical reaction, the system tends to 
approach chemical equilibrium where the forward and reverse reaction rates are the 
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' ► 

Extent of Reaction 

Figure 3. Non-ideal Mixture: G vs. ^ 

same and concentrations of the reactant and product species do not change with time. 
When the equihbrium condition is achieved, proportions of the various compounds 
remain unchanged, and the reaction ceases to progress. Prior to reaching the point of 
equihbrium, the system fluctuates between different equihbrium microstates. Suppose 
that the system is confined in a lower (more stable) Gibbs free-energy minimum and, 
occasionally, a fluctuation may be large enough to push the system over the maximum 
to the region of higher energy, i.e. a matastable minimum. A small fluctuation can 
overcome the shallow barrier back to the more stable equilibrium state [4]. Any 
thermodynamic system, in this case a chemical reaction, tends to eventually reach 
the lowest minimum in the Gibbs free energy. Naturally, if the "unstable" barrier 
is too high or the minima are far apart a shift of the equilibrium from one local 
minimum state to another is less probable. 

Within this picture, the local curvature of the Gibbs free energy is positive for 
all points except those between the two inflection points. Moreover, the portion of 
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the curve between the minima at tlie inflection points is said to be locally stable but 
globally unstable. In this region on the curve metastable states occur which appear 
to be stable to small perturbations, but mixed configurations at the same extent of 
reaction represent more stable states with lower free energy. A straight line connecting 
the two minima corresponds to a phase boundary [4], i.e. a phase transition from the 
phase at one minimum to the phase at the other minimum. Positive local curvature 
fails at the points of infiection. Local stability determines whether, after a small 
perturbation, a system will return to the original equilibrium state. Here our focus 
is on examining conditions leading to failure of local stability. 

Remcirk 4. Ideal Mixture. For an isothermal and isobaric ideal mixture, 

dArG _ RT 

When T 7^ 0, consider that 

which is the so-called logistic equation [2]. The process described by this equation 
has two equilibrium positions, namely ^ = and ^ = 1. Between these two points the 
field is directed from to 1. As a result the equifibrium position ^ = is unstable (as 
soon as the reaction proceeds away from ^ = reactants are converted to products). 
Meanwhile the equilibrium position ^ = 1 is stable. Moreover, integral curves tend 
asymptotically to the line ^ = 1 as Aj-G — > +00 and to the line ^ = as A^G —00. 
Such curves describe the passage from one state (0) to another (1) in an infinite A^G. 

3.2. Compounds. Consider a generic chemical reaction written as, 

Wb\B + ... |^'5|'S'+ \ut\T + ... 

where components on the left side are designated reactants while components on the 
right hand side are products [9]. The lu^l, i — A,B,... are stoichiometric coefficients 



(3.20) 
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of the reaction. Another formal representation of a chemical reaction which better 
lends itself to mathematical manipulation is given by [9] , 

lysS + utT + ... + vaA + ubB^J2 = (3-22) 

i 

The chemical potential of each component is given by [9] , 

IJLi = li*{T,p) + RT\Yi-iiXi = iJL*{T,p) + RThiai i = A,B,... (3.23) 

where a, is the activity of component i. Recall that Xi = W^^, with — + 
and that 7^ depends on the extent of reaction. Thus, the Gibbs free-energy of the 
reaction is given by [9], 

I^rG^Y.ii.v^ = Y, v^ix\ + it!T In fj = + In J] 7^4" (3.24) 

i i i i i 

Introducing the quotient of reaction, Qa — Jli '^T ~ Y\i iT Y\i ~ QiQci where the 
subscript c denotes concentration, the expression (3.24) can be written as [9], 

A^G = A^G^ + RT In g„ (3.25) 



where A^G^ = i'iH*{T,p^) is the standard Gibbs free energy of reaction. It follows 
that, 

QaiO = (3.26) 

Naturally, if the system reaches equilibrium, namely Aj.G = 0, the parameter Qa 
is denoted by K^i the equilibrium constant, and the expression in (3.26) becomes [9], 

= X„ = (3.27) 



Thus, for an isothermal and isobaric single chemical reaction. 
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^^E(|) - + 1 1^(11 ^D] (3.28) 



where 



For simplicity, denote ^ln(]J^7]'') = Then, the expression in (3.28) can be 

rewritten as, 

^ = -«^E(|)-^|^ + WK)l (3-30) 

The expression in (3.29) denotes the influence of the relative amounts of reactants 
and products at each extent of the reaction while W{^) represents the relative strength 
of non-ideal (inter-particle) interactions existent between products and reactants. Con- 
sequently, at any value of the extent of reaction, there is a "possible" value of W such 
that the two mentioned forces exactly balance one another. Thus, at a given certain 
extent of the reaction determined by the relative amounts of reactants and products, 
W at that point corresponds to the relative strength of non-ideal interactions that 
must exist between products and reactants for a failure of local stability. 

Remark 5. For the ideal gas mixture, Qa — Qc— i4| and therefore, ^ = -jjz^ > 
(see Fig. 4). In this case, — 1 and the expression in (3.30) reduces to, 

f -«^lE(|)-§fl (3.31) 

Moreover, when the sum of the stoichiometric coefficients vanishes (i.e the isothermal- 
isobaric Ideal gas mixture, see (3.19)), the expression in (3.30) reduces further to. 



(3.32) 
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0.2 0.4 O.B 

Quotient of reaction vs. eKtent of reaction - Ideal MiKture 



Figure 4. Ideal Mixture 



Theorem 2. Let T ^ 0. Then, for an isobaric and isothermal single chemical reac- 
tion, = if and only if 



w{0 



(3.33) 
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The Gibbs free energy, G, is a convex function of the extent of reaction whenever, 



The curve described by the expression in (3.33) is denoted as the curve of phase 
boundary in the W-^ plane. Such a curve traces the phase boundary between the 
convex and concave regions of the Gibbs free energy. In particular, at a fixed value of 
the extent of reaction, the system is locally stable whenever the value of W is less (in 
absolute value) than the value on the curve of phase boundary. If instead the value 
of W is greater, the system is locally unstable. 

In early stages of the reaction, reactant species are far in excess of product species, 
and W is (in absolute value) relatively large. Since a change from reactant phase 
to product phase is improbable early in the reaction, interactions between products, 
favoring product formation, must be greater than those between reactants, favoring 
the reactant phase. W indicates the balance of the strengths of the product and 
reactant interactions required for failure of local stability. As the reaction progresses 
toward the extremum of the curve of phase boundary, ^ = 0, the relative difference 
in strength between the two types of interactions is a minimum. At this critical 
extent of reaction, a change between reactant and product phases requires the smallest 
difference between their constituent interactions and is thus most probable. Past this 
critical point the extent of reaction increases. To achieve local instability the relative 
strength of the interactions favoring reactants must be increasingly greater than those 
favoring the products. 




(3.34) 



and a concave function whenever. 




(3.35) 



Example II: Synthesis Reaction 
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Consider a simple synthesis reaction in which two or more substances combine to 
form a more complex substance. For example, 2 moles of di-hydrogen react with 1 
mole of oxygen to give 2 moles of water, 

2H2 + 02^ 2if20 (3.36) 

with 

7V^, = 2-2^ No,^l-i and Nh,o ^ 2^ 
and the corresponding stoichiometric coefficients given by, 

uh^ — —2 1^02 — ~1 ciiT'd ^HiO — 2 

Then, from (3.33), 

6 

^^^^ ^ "e(l- 0(3-0 

A plot of the curve of phase boundary for this reaction is given in Fig.5. 

The extremum of this curve defines the "most probable" transition point between 
the reactant and product phases. At this point differences in the relative amounts of 
reactant and product species, and relative differences in the strengths of their non- 
ideal interactions, are minimal. Taking the derivative with respect to ^ (Fig. 6) and 
setting it to zero yields, 

i = 0.4514 

It follows that W — —9.5. The fact that the critical extent of the reaction is less 
than 0.5 suggests that the products and associated non-ideal interactions are more 
strongly favored such that the product phase is preferred even before half the extent 
of reaction is reached. 



(3.37) 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

extent of reaction 



Figure 5. Synthesis Reaction: W vs. ^ 

Now in analogy, consider the dissociation reaction, i.e. the synthesis reaction in 
the opposite direction. In this case 2 moles of water split into 2 moles of hydrogen 
and 1 mole of oxygen. Namely, 



2H2O 2H2 + O2 



(3.38) 
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Figure 6. Synthesis Reaction: W vs. ^ 



Following the same steps as for analysis of the synthesis reaction, the following 
expression for W is obtained, 
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0.2 0.4 O.B 0,0 

extent of reaction 



Figure 7. Dissociation Reaction: W vs. ^ 

The graph of this curve of phase boundary is shown in Fig. 7. Taking the derivative 
with respect to ^ and setting it to zero (see Fig. 8), 



^ = 0.5486 
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0,3 0,4 0.5 

extent of reaction 



Figure 8. Dissociation Reaction: W vs. { 



Note, W = —9.5 is apparently "invariant" to the direction of the chemical reaction. 
The critical extent of reaction, ^ = 0.5486 for the dissociation reaction, consistent with 
inferences drawn from analysis of the synthesis reaction. That is, interactions between 
the synthetic species {H2O) are more favorable than those between the individual 
species O2). The sum of the two critical extents of reactions is 0.4514 + 0.5486 = 
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1. 

Example III: Single Displacement Reaction 

A single displacement reaction is one in which an atom (or ion) of a single com- 
pound replaces an atom of another compound. As an example, consider the single 
displacement in which copper ions in a copper sulfate solution are displaced by zinc, 
forming zinc sulfate: 

Zn + CuSO^ ^Cu^ ZnS04 (3.40) 

with, 

Nzn^i-^ NcuSO,^i-C Ncu^i and Nznso.^i 
The corresponding stoichiometric coefficients are given by, 

i^Zn = -1 l^CuSOi = -1 t'Cu = 1 and PZnSOi = 1 

Then, the expression in (3.33) is given by. 

The graph of this curve of phase boundary for this reaction is shown in Fig.9. Taking 
the derivative with respect to ^ and setting it to zero (see Fig. 10), 

e = 0.5 

It follows that 1^ = — 8. In this case ^ = 0.5 corresponds to the minimum difference 
between interactions associated with the product and reactant species for a phase 
transition to occur. Note that, for the displacement reaction, the critical extent of 
the reaction is 0.5. In this case, the number of products species equals the number of 
reactants species and the critical extent of the reaction is independent of the direction 
of the process. This is in contrast to what was found for obtain in the synthesis and 




0,1 0.2 0.3 0.4 0.5 O.B 0,7 0.0 0.9 



extent of reaction 



Figure 9. Single Displacement Reaction: W vs. ^ 

dissociation reactions where product and reactant species are not equal and the critical 
extent of the reaction depends on the reaction direction. 

3.3. Multi-chemical reaction. Here the Gibbs metric is introduced for a closed 
system with r chemical species in which / independent reactions can occur. The total 
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Figure 10. Single Displacement Reaction: W vs. ^ 

change of mass, drrii, is equal to the sum of the changes resulting from the different 
reactions, [16], 

drrii = Mi^Uind^n n = l,..,l (3.42) 

n 
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The principle of Conservation of Mass can be stated as [16], 

(^^^J2Y1 ""inMidin = (3.43) 
i n 

As done previously, consider the number of moles of the components in the system 
instead of their masses. The change in number of moles of component i is given by, 

dNi = ^ l^ind^n 
n 

The Legendre submanifold Sq is defined by the constitutive relation G — G{T, p, ^i, ... 
and the expression in Eqn.(3.3) becomes, 

dNi,n ^ ^ dNi^ ^ 

where ^„ is the extent of the n-th reaction with i 
consider the differential of the Gibbs free-energy. 



... = d^n 

— 1, ...,r and n — 1, I, [9]. Now, 



dG = -SdT + VdP - A^din 



(3.44) 



By recalling the expression in (1.16), the following metric is obtained. 



T 



Vijc 



aV 

aV -krV 

ArSi ArVi 



-ArSi 



ArVr 
'dAi\ 

dA, 



(3.45) 



This is the Gibbs metric for a multicomponent thermodynamic system in which I 
independent chemical reactions occur. 
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4. System 3: Open systems 



Finally, consider open systems in the absence of external fields. Recall that the con- 
tact 1-form 9, restricted to the Legendre submanifold Sg described by the constitutive 
relation G — G{T,p,Ni, ...,Nr), gives the following differential, see (1.39), 



dC = -SdT +Vdp + J2 l^idNi 

i 

Define a partial molar quantity as [9] , 



(4.1) 



(4.2) 



Since the Gibbs metric is defined as the Hessian of the Gibbs potential, rjij^ — g^ig^j , 
where and are the extensive variables, the following result is obtained. 



Vijc 



T 

aV 
-'Si 



aV -Si 
-krV Vl 



\^ — Sr Vr A*' 



rl 



-sl\ 

Hrr J 



(4.3) 



where //^fc = {§^)t,p,n^^,- 

At constant temperature and pressure, conservation of energy (4.1) reduces to. 



dC = ^ HidNi 



(4.4) 



Therefore, the metric in Eqn.(4.3) becomes. 



( A*ii A«2i 



(4.5) 
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4.1. Ideal solutions. Now, introduce tlic Gibbs metric for tlie cases of an ideal and 
non-ideal solution in an open thermodynamic system without chemical reactions. A 
component i in solution is said to be ideal when its chemical potential is given by, [9], 

lii^ fil{T,p) + RT In Xi (4.6) 



where ii*{T,p) is the standard state chemical potential which is independent of the 
composition and Xi — = ^ is the mole fraction. Note, the total number 



of moles, N, depends on the number of moles of the r components, namely N — 
N{N,,...,Nr). 

Using the expression (4.6), the metric for an ideal solution that depends on the 
temperature and the pressure of the system is obtained. 



Vijc 



( 



aV 



+Rlnxi 



dT 



v. 



Ni 



aV 
-krV 



,-9^)T,Ni 



[^)T,Ni 



RT 
' N 



'rwr)p,N, + i?lnxA 



^ dp ' 



RT 
' N 



(4.7) 



For an isothermal-isobaric system such a metric becomes. 
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(4.8) 



This metric is degenerate everywhere since det(77yQ) = for all values of the number 
of moles of the r components. This result indicates the Gibbs free energy metric for 
an open system does not distinguish differences between ideal components. As far as 
the metric is concerned, due to the explicit lack of inter-component interactions, the 
system is apparently comprised of a single, indistinguishable, ideal component. 
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4.2. Non-Ideal solutions. For non-ideal systems, the chemical potential of compo- 
nent i takes the form, [9], 



^ii = l^*{T,p) + RTln'yiXi (4.9) 



where 7^ is the activity coefficient of species i that depends on temperature, pressure 
and composition of the solution. Expression (4.9) considers deviations from ideality. 
The partial molar entropy. Si, and the partial molar volume of component i, Vi, are 
given by. 



and 



Moreover, 



and. 



-S. = + Rl\m,Xi + T(^)„„.] (4.10) 



K^(f)..+«T(^),„, (4.11) 



^^ii - {-^)T,p,N^^i - ~ TV-* ( Q^')T,p,N^-t] (4.12) 



, 1 dlnji 



46 



MANUEL SANTORO AND ALBERT S. BENIGHT 



For symplicity, denote {^^^)T,p,Nj^k ~ i^^'yi)k retain the usual notation for 
partial molar entropy and volume. Substituting the expressions (4.10)-(4.13) into the 
metric (4.3), the general metric of a multicomponent non-ideal solution is obtained, 
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(4.14) 



and, in the case of an isobaric-isothermal non-ideal system, the Gibbs metric simplifies 
to. 



/[(ir-^) + (ln7i)i] [-^ + (ln72)i] 



[-^ + (ln7^i] \ 



V [-^ + an7.)i] [-^ + (ln7,),] 



[(i-^) + an7.W 
(4.15) 



= RT[ 



' N 



({— - -) 

AT VaTz N> 



\ 



'at 



' N 



N 



+ 



/^(ln7i)i (11172)1 ■■■ (ln7r)i\ 
(In 72)1 (In 72)2 ••• (ln7r); 



(i-^)/ \an7r)i (ln7r-)2 



(ln7r),7 
(4.16) 



{.Vija^ideal ~l~ {jlijc) deviation 



(4.17) 



where the subscript " deviation" indicates the degree of non- ideality. 
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5. Conclusions 

Geometrical thermodynamics has been apphed to study the behavior of several 
simple thermodynamic systems. For a single component, closed system intriguing 
relationships between the geometrical concepts of degeneracy and scalar curvature 
of the Weinhold metric, and the physical concepts of a phase transition and inter- 
particle, non-ideal interactions are divulged. For multi-component closed systems the 
Gibbs metric was presented and the scalar curvature and degeneracy of this metric 
was determined for a few simple cases and could be indicative of physical behavior. 
In summary, this study provides convincing examples that this geometrical approach 
to analysis of thermodynamic systems can be applied to actually divulge important 
physical and chemical behavior. 
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